Can gravitation accelerate neutrinos? 
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I. INTRODUCTION 

The well known September 2011 OPERA neutrino ex- 
periment [l| has produced a myriad of articles trying to 
explain or understand, prove or disprove the experimen- 
tal results for superluminal neutrino motion. One of the 
most important contributions to the understanding of the 
problem is the Cohen-Glashow work [2|- 

Most of articles are based on models which, by con- 
struction, ignore or neglect general relativistic effects. 
Furthermore, even though the general relativistic aspects 
of the dynamics are not considered, the arguments pre- 
sented in these papers seem to tacitly assume that neu- 
trinos follow geodesies in the presence of gravitational 
fields. 

It has been known for quite some time [3|-l6[ that spin- 
ning massive test particles (tops) follow non geodesic 
paths when moving on gravitational fields. 

Even though a recent communication seems to refute 
the fact that the original OPERA experimental results do 
indeed describe superluminal propagation of neutrinos, 
we believe that is interesting to show that a consistent 
general relativistic model allows for the possibility of su- 
perluminal propajgation of tops moving on gravitational 
backgrounds [a-0]- 

In this article we discuss an exact solution for the equa- 
tions describing the motion of a top on a Schwarzschild 
background. A Lagrangian derivation for these non 
geodesic equations of motion obeyed by tops moving on 
agravitational background was first obtained by Hojman 
@, m (using a formalism by Hanson and Regge Q devel- 
oped for flat spacetime) which it is outlined in the next 
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section. We show that massive tops may reach super- 
luminal velocities provided that their particle mass m, 
energy E and total angular momentum magnitude j of 
the orbit, satisfy a relation (which is explicitly derived in 
Section IV). 

Non geodesic equations of motion for tops were first 
derived by Mathisson Q and Papapetrou Q as limit- 
ing cases of rotating fluids moving in gravitational flelds. 
On the other hand, massless spinning particles do follow 
null geodesies as showed by Mashhoon Q who used the 
Mathisson-Papapetrou formalism for his derivation. 

Loosely speaking, in the presence of gravitational 
fields, massive spinning particles (such as neutrinos) fol- 
low non geodesic paths [sHa] while photons move in null 
geodesies (in spite of their spinning nature) because they 
are massless 0. 

Even though the Equivalence Principle is sometimes 
interpreted as stating that test particles in a gravita- 
tional field should follow geodesies, this interpretation is, 
strictly speaking, valid only for spinless point test parti- 
cles. Extended particles arc, in general, subject to tidal 
forces and follow, therefore, non geodesic paths. 

In the Lagrangian formulation of the motion of tops, 
the velocity u^ and the canonical momentum P^ vec- 
tors are, in general, not parallel. For the motion of tops 
in the presence of electromagnetic and/or gravitational 
fields, the (square of the) mass m?{= P'^P^ > 0) is con- 
served implying that the momentum vector remains time- 
like along the motion. Nevertheless, the velocity vector 
may become spacelike [3, [^ [ill [l3l ■ 

The proper treatment of the (subtle but crucial) lack 
of parallelism between velocity and momentum is best 
achieved with a Lagrangian formulation of the motion of 
tops, because otherwise the canonical momentum cannot 
be appropriately defined. Furthermore, the Mathisson- 
Papapetrou formulation gives rise to third order equa- 
tions of motion, while the Lagrangian approach gives rise 



to second order ones [3, [Tol ■ 

It turns out that for suitable choices of initial condi- 
tions this non geodesic motion can accelerate tops (and 
neutrinos) beyond the speed of light. 

Theoretical results involving superluminal propaga- 
tion of massive spinning particles and fields in inter- 
action with electromagnetic or gravitational fields have 
been reported previously in the literature by Velo and 
Zwanziger [11|, Hanson and Regge Q, Ho j man [5| and 
Hojman and Regge jl2| . for instance, while experi- 
ments reporting hints of superluminal neutrino prop- 
agation can be found in Refs. [ij, ^^, among oth- 
ers, and other results reported in the MIT webpage on 
Superluminal Neutrinos (The Net Advanc e of Physics, 
http://web.mit.edu/redingtn/www/netadv I . 



II. MOTION OF A RELATIVISTIC TOP ON A 
GRAVITATIONAL FIELD 

Consider a relativistic (spherical) top. Denote its po- 
sition by a four vector x^ while its orientation is defined 
by an orthonormal tetrad e(c)^- A gravitational field is 
described as usual in terms of the metric field 5^^ 0, Q . 
The tetrad vectors satisfy 



with 77(„^)(= Ty^"^)) given by 

Viad) = diag (4-1,-1,-1,-1) 



(1) 



(2) 



and have, therefore, six independent components (con- 
sistent with the number of parameters of the Lorentz 
group). 

The velocity vector u^ is defined in terms of an arbi- 
trary parameter A by 



dxi" 



(3) 



Besides, the antisymmetric angular velocity tensor cr^'^ 



<J -V e(„) ^^ 



1//^ 



(4) 



where the (covariant) derivative De^i^)'^ /DX is defined in 
terms of the Christoffel symbols V^ pr , as usual, by 



De 



(/3) 



de 



W) 



DX 



dX 



+ ^''pT e(i3f u 



P iJ 



(5) 



Note that general covariance may be achieved unam- 
biguously at the level of the Lagrangian formulation [5| 
because only first derivatives of the dynamical variables 
are used in the construction of the Lagrangian. If im- 
plemented using the equations of motion, terms propor- 
tional to the Riemann tensor may be missed, because 
(flat spacetime) second partial derivatives commute while 
(curved spacetime) second order covariant derivatives do 



not commute and their commutator is proportional to 
the Riemann tensor. If no Lagrangian theory for a sys- 
tem of special relativistic equations of motion is known, 
the introduction of gravitational interactions cannot be 
unambiguously implemented. 

The Lagrangian L = 27(01,02,03,04) is constructed 
as an arbitrary function of four invariants 01,02,03,04 
such that the action S = J L dX, be A-reparametrization 
invariant 

L(ai, 02, 03, 04) = (oi)^/^/: (02/ai, 03/(01)^, 04/(01)^) , 

(6) 
(the speed of light c is set equal to 1, at this time), C is an 
arbitrary function of three variables and oi = u^^Up, 02 = 
<y^'^<y^iv = — tr(f^), 03 = Uaa^^ap^u"^, 04 = dct{a). 
Note that this expression seems to require oi to be pos- 
itive. To be precise, the fact that oi be positive is not 
crucial in the formulation, as it will be proved later on 
[see Eq. (|2ip ]. One could rewrite the expression for the 
Lagrangian as 

L(oi,02,03,04) = (02)^/^^1(01/02,03/(02)^,04/(02)^), 

(7) 
for instance. 

The conjugated momentum vector P^ and antisym- 
metric spin tensor S^i, arc defined by 



dL 



^au 



dL 



^h-fi ■ 



(8) 



(9) 



The equations of motion are obtained by considering 
the variation of the action S with respect to (ten) inde- 
pendent variations Sx^ and (the covariant generalization 
of) Se^"' defined by 



^"=r;("'^)e(„)''(5e(^)'' =-(50"'^ 



(10) 



The (non geodesic) equations of motion turn out to be 

m 



DPP- ^ 

DX " 2 ''"^ ' 



(11) 



and 



DX 



= S^'^ax" - cj^^Sx" = P^u" - u^'P'' . (12) 



These results hold for arbitrary C. The dynamical vari- 
ables P'' and S'^'^ may be interpreted as the ten gen- 
erators of the Poincarc group. In order to restrict the 
spin tensor to generate rotations only, the Tulczyjew con- 
straint [la 



S^"'P^ = . 



(13) 



is usually imposed |a, Q- R turns out that both the 
top mass m and its spin J are conserved quantities (see 
Appendix A) 



m^ = PI'P,, 



(14) 



72 — _ qfj.1^ q 



Furthermore, if ^'' is a Killing vector, then 



1 



Q EE P^Cm - -S^-'C^,. ., 



(15) 



(16) 



is a constant of motion 0, [a, [T^l ■ 

The following gauge choices and "invariant relations" 
(defined in Appendix B) 



c° = A = t . 



e^of^P'^/m, (17) 



may be implemented (5|, |6|, |8| to fix the arbitrary param- 
eter A and to restrict the (Lorentz transformations) six 
degrees of freedom of the tetrad to three dimensional ro- 
tations (for details, please see Appendix B). The previous 
choices satisfy condition P^ 



P^Pu 



e{^,^'^e|m,_ = tA = ^7(00) = 1 



'(0) e(o)^ 

and arc consistent with constraint (|T3 
DiS^'^P,) 



(18) 



DA 



= 



= (S^Va"^ - a'^^^A'') P. 



lS^"'- 



De 



(0)^ 



DX 



- mS'^^ax^e^o), - mS'^" G^xei^^)^ - (19) 
where we have used the fact that 



DX 



-<^^\e-{a) , 



(20) 



for any of the tetrad vectors, according to the Eqs. (|H), 
(P) and ©. 

The consistency of the constraint (|13p with the equa- 
tions of motion (|lip and (|12p is guaranteed by making 
use of the arbitrariness of Lagrangian C (or £1) in ^ 
(or (O) by appropriately constructing it. As a matter of 
fact, the Lagrangian Q 



L 



Aai — Ba2 



-^/{Aal - Ba2Y - 8B{Aa3 - 2^04) (21) 



1/2 



gives rise to the equations of motion ([TT|) and ([T^ and the 
Tulczyjew constraint (fT3| plus a Regge trajectory defined 
by Bm'^ — AJ'^/2 = AB. Besides, it can be proved this 
Lagrangian is well defined for ai < 0. Therefore, one 
may consider Lagrangian ([21]) as the starting point of 
this theory. 

In order to prove that ((2T|) is a well defined Lagrangian 
for any value of Oi, define 



D = -^/{M^[^^B^i:2fl^T6B^^M^^MB^, (22) 



such that the Lagrangian ((2T|) now becomes 



(23) 



In what follows we show that D > {Aai — Ba2)/2 which 
ensures that L is real for any value of ai under very 
general assumptions. 

First, we can see that 04 is always non negative because 
the determinant of an antisymmetric matrix in an even 
dimensional space is a perfect square. The sign of 03 can 
be found as follows. We have that 






may be rewritten as 



0-3 



Uo^U", 



ith 



TjOi _Q!/3 

U = a ' u^ , 



(24) 



(25) 



(26) 



Consider the momentum vector P" which is always 
timclike (due to mass conservation). We now show that 
Ua is a spacelike vector because it is orthogonal to P". 
The orthogonality relation reads 



DPf^ 



(27) 



where we have used ([201) for e^^x = P'^/m, and the fact 
that Up is orthogonal to DP^ /DX as it can be easily seen 
from ([TT|) . Therefore, UaU"" < 0, and then 03 is positive. 
It is enough to choose AB < to end the proof that the 
Lagrangian ([^T]) is well defined irrespective of the sign of 
ai. Even if AB > there are regions where Lagrangian 
([^ is wcU defined for oi < 0, but for AB < 0, it is well 
defined everywhere. 



III. EXACT SOLUTION 

The equatorial motion of a top in a Schwarzschild field 
background may be solved exactly. The Schwarzschild 
line clement is ds^ = gttdt^ + grrdr^ + geedO^ + g^^dcfP' , 
where gtt = c^{l - 2ro/r), g„ = (1 -2ro/r)"\ ggg = 
^^^, g<t><t> = — ?'^sin^6', and 2ro is the Schwarzschild ra- 
dius. From now on we reinsert explicitly c in all the 
expressions. The general equations (jlip . (|12p were writ- 
ten in Ref. H along with ([T^ . ([Ti j) . ^. ([TC j) and ([T7 |) 
for the four Killing vectors of the Schwarzschild metric. 
In this article, we restrict ourselves to the motion in the 
plane defined by cos = 0. If the top is initially in that 
plane and = 0, then it remains in the equatorial plane 
[3, in which 9 ~ 7r/2 and P^ = 0. This reduction is 
possible due to the fact that the direction of angular mo- 
mentum is conserved, so two of the four Killing vector 
conservation laws are used to restrict the motion to the 
equatorial plane. 



It is convenient to define the dimensionless parameter 



jVo 



(28) 



where J = h/2 is the top's spin (as well as the neutrino's 
spin). Thus, the set of equations ((TT|) - ((T7l) (including the 
two remaining Killing vector conservation laws, energy E 
and total angular momentum magnitude j in addition to 
the conservation of mass and spin) may be solved exactly 
to yield [| 



Pt 



-j±EJ/{mc^ 
1 — 77 



ETJJro/{mr^) 
1-r, 



(29) 



(30) 



and, from Pf^P^ = rri^c^, we get 



P'' = ± 



P' 



Pi 






1- 



2ro 



-,1/2 



(31) 



We can now solve for the velocities. To this end, we 
use two of the equations of motion p^ . In the plane 
defined by = 7r/2, these equations become 



DX 



P*r - P'' 

S'f"'PA, DPt DS''"' Pa, 



and 



DX 



Pi DX 



= P'cf.-P'^ 

S'f"'Pr DPt 



DX Pt 



S't"' DP^ 
Pt DX 



.(32) 



DS't"' Pr S*"" DPr 



P? 



DX 



DX Pt Pt DX 



,(33) 



which may be solved for f and 0. To perform this task, we 
use the equations of motion ([TT|) and ([T^ . the relations 
between the spin and momentum [5| 



S'' 



gti 



S'>"'P^ 

Pt 



S'f"'Pr 
Pt 



(34) 



(35) 



which are consequences of the constraint ((T^ . and the 
condition 



(S"^"^) 



2 .P {Ptf 



(36) 



which is derived from Eqs. ([T51) . ([Ti)) and P^ . 



Taking these results into account, the velocities turn 
out to be 



^^'> 



2ro\ f 2r^ + l \ f P^ 

v-ij \Pt 



2ro 
r 



Finally, we get 



dr 



277 + 1 



1] 



1 



pr 

'P't 



Pi> 

,,2pr 



(37) 



(38) 



(39) 



It is worth noting that the three preceding expressions 
coincide with the usual results for geodesic motion when 
J^ = (and therefore 77 = 0). Once the solutions are 
spelled out, we can find one of the main results of this 
work. From expressions ([57)1 and ([55)) . we find 




777 



n2(l-A) 



(P*) 



(40) 

where c^P* = (1 — 2ro/r) ^Pt, and we define the super- 
luminal parameter A as 



A 



3r/(2 + 7;) 
m'^c^r^^l — 77)' 



:iP4.f>0. 



(41) 



From (j40|) is straightforward to realize that it is possible 
that ds^ < for some of the solutions, at least in part 
of the top's trajectories. It is clear that the contribution 
of the A parameter is important for small mass parti- 
cles, such as neutrinos, because the A dependence on the 
particle mass behaves as 777"'*. 

In the next section we study the trajectory of small 
mass particles under the approximation in which the top 
moves far away from the Schwarzschild horizon, i.e., for 
r ^ rg. We find the conditions that the mass, energy 
and total angular momentum of the top must satisfy in 
order to produce superluminal motion in a segment of its 
trajectory. 



IV. SUPERLUMINAL MOTION FOR r > ro 

The main purpose of this section is to show that there 
are some particle trajectories such that A < 1 in part 
of the path and A > 1 in the rest of it. Knowing that 
J = h/2 is the spin of the top, superluminal behavior 
imposes a condition on the mass 777, the total energy E, 
and the total angular momentum j. 

We look for a solution such that the particle is mov- 
ing at distances r much larger than the Schwarzschild 
radius of the black hole, rg/r ^ 1, i.e., we only consider 
the motion of the particle in a weak gravitational field. 
Besides, to focus in regimes in which A is relevant, we re- 
strict ourselves to the case of a small mass particle, such 



that its total energy is larger than its rest mass energy, 
E ^ mc^. Both assumptions imply that ij <ti 1. 

In the weak field and small mass approxima- 
tions, taking the upper signs in the preceding so- 
lutions p9|) and (pO)) . the momenta become P^ = 
{l + ri){-J + EJ/mc^), Pt = il + T]){E-jJro/mr^), 
and 



P'- = ± 



E' 



1 



-^-z^\J 



EJ 



-,1/2 



(42) 



In order to study the particle motion, we assume that 
initially the particle approaches the central body from 
infinity and it remains always at a distance r ^ tq. Thus, 
the orbit of the particle in the equatorial plane can be 
characterized by two values of r in the trajectory, r^ and 
r^. The first one is the value of the orbit's radius in 
which the particle's velocity reaches the speed of light, 
i.e., the point in which A = 1. The second one is the 
return point of the particle's orbit, and it is defined as 
the point when its radial momentum vanishes, P^ = 0. 

By imposing the condition A = 1, it is straightforward 
to obtain that the value of Tc is given by 



6JVo 



EJ 



mc^ 



1/5 



(43) 



On the other hand, we can calculate the return point 
for the particle's trajectory r^ solving /"" = from (|^. 
We get 



V. CONCLUSIONS 

We present an exact solution to the non geodesic equa- 
tions of motion in a Schwarzschild gravitational back- 
ground that allows for superluminal propagation of mas- 
sive spinning test particles. The superluminal motion de- 
pends strongly of the inverse of the mass of the particle. 
This is, of course, relevant for small mass particles such 
as neutrinos. We have shown that a consistent general 
relativistic theory which allows for superluminal propa- 
gation is possible. 

Furthermore, the superluminal propagation effect pre- 
sented here can be achieved in weak gravitational fields, 
as for example, in the surface of the Earth. In the pres- 
ence of stronger gravitational fields (which can be easily 
found in astrophysical context) this effect will be much 
enhanced. 

It is worthwhile mentioning that the aforementioned 
effect depends strongly on the fact that we deal with 
small mass particles, as the relevant A parameter has an 
TO"^ dependence, and it is therefore very unlikely that 
superluminality could be detected for particles other than 
neutrinos. 

Finally, we would like to mention that the equation 
of motion (|lip can be rigorously generalized to include 
the gravitational self- force of the tops [la] , showing that 
these corrections also modify the geodesic paths. How- 
ever, the new forces are proportional to the particle mass 
and therefore their effects for small mass particles are 
negligible compared with those presented in this work. 



rR 



21 
E 



1 



JE 
j mc^ 



(44) 



The particle reaches the speed of light in re (be- 
cause A(rc) = 1 and therefore ds^ = 0). However, to 
have superluminal motion after that point requires that 
A(rfl) > 1 (implying ds^ < 0). Using Eq. ^ in (gT]), 
this implies a condition that the mass, the energy and 
the total angular momentum of the particle must fulfill 



6J^^Vo > m'^c^ 



J 



EJ 



(45) 



The above expression shows the relation that m, E and 
j satisfy to achieve superluminal motion. Notice that the 
condition cannot be met by spinless particles. Knowing 
the energy of the particle, it is possible to use (^5]) to 
estimate the mass and total angular momentum of the 
particle if superluminal motion is ever detected. 

Interestingly, one can prove that condition (j45p is 
equivalent to tr < re- This means that when the particle 
has a "ballistic" trajectory around the black hole, it can 
always be superluminal in some part of its orbit, nearest 
to 2ro, for appropriate m, E and j. Clearly its speed will 
not be constant along the trajectory, because in the first 
part of the trip, the gravitational field speeds it up until 
it reaches r/j, slowing it down afterwards. 



Appendix A 

The mass conservation law may be obtained as follows. 
Rewrite the velocity vector u^ as 

u"" = -^{P^PnW' + ^iP^u^^-P^unP'' 

= i-(pA'u--p-wA<)p^ + i-p^u'^p-. (Al) 

We use the equation of motion for the spin tensor (fT2)) 
to deal with the first term on the right hand side of the 
previous equation. Then 



m^u" = P, 



DS^"' 
DX 



P^.u^^P" , 



DP 

-—^S'^-' + P^^u^P''., 



(A2) 



where we have used the constraint S^^P^ = 0. We multi- 
ply by DxPv, and due to the fact that DxP^,S^"' DxP^ = 
0, we find 



DP 

m u — = PnirP 



.DP, _ 



DX 



(A3) 



However, using the equation of motion for P,, we find 
that u'^D\P, = 0, and then we get the condition 

P-^ = 0, (A4) 



which imphcs that m^ is constant. 



Appendix B 

In this Appendix we show how to implement the 
choices 

e^of = P''/m. (Bl) 

using a different approach to the one presented in [8[. 
Let's start by defining 8(^^13) by 

S{ap) = e(a)^S'^''e(^)^ . (B2) 

It is a straightforward matter to reahze that the six quan- 
tities Stap) s-rc constants of motion 



DS. 



(q/3) 



DA 



= 0, 



because of (fT2l) and (l20l). We can now choose 



S, 



m 



0, 



(B3) 



(B4) 



as three initial conditions which are, of course, preserved 
in time because of (|B3[) . Hanson and Regge call these 
conditions "invariant relations" . 

Conditions (|B4p imply that the vector e(o)/i is a null 
eigenvector of the spin matrix S^'^ , i.e.. 



S'^-'t 



(0)^ 



= 0. 



(B5) 



The spin matrix S^'^ has even rank (because it is anti- 
symmetric). Due to the fact that S^'^Pi, = 0, its rank is 
not four, so it must be two (otherwise it would be zero, 
rendering it trivial). If its rank is two, it must have two 
null eigenvectors, which are the momentum vector P^ 
and the Pauh Lubanski vector W = ^e^'''°''^ SapPy, as 
one can easily prove (see also Q)- 

Therefore, the vector e(o)'' may be expressed as a linear 
combination of the momentum vector P'^ and the Pauli 
Lubanski vector W^^ 



PM 



Hoy 



TW^' . 



(B6) 



Hanson and Regge [8| construct the Hamiltonian the- 
ory of the top (pages 523 and following of reference Q) 
using Dirac's method [l3l- There they handle the Tul- 
czyjew constraint S^'^Py = by extracting its first class 
content $2 (in Dirac's terminology) with 



cf>2 = ^e^''"''V^«/3«0, 



(B7) 



where the sign "« 0" is read "weakly equal to zero". 
This means that the constraint $2 vanishes, but its Pois- 
son bracket relations with some dynamical variables are 
different from zero. Loosely speaking, Dirac's method 
provides one way to mend this contradiction by redefin- 
ing the Poisson brackets relations. The new brackets are 
called Dirac brackets in his honor. The Dirac brackets of 
the constraints (with any dynamical variable) arc identi- 
cally zero. 

Due to the fact that $2 is first class, it generates gauge 
transformations [l7|. Therefore a gauge (associated to 
it) may be chosen. If the Poisson bracket of a dynamical 
variable A with $2 is such that [A, $2] ~ 0, then A is 
invariant under $2, i-c, is gauge invariant. So, in order 
to choose a gauge associated to $2 one needs to find 
a variable B such that its Poisson bracket with $2 be 
different from zero. 

Use dm to get p 



3(0) 



P 



(B8) 



and consider the Poisson bracket relations, equations 
(3.11) of reference Q 



\p, -.f^ Q'^P] — p, "nM0 _ p, ./SnA'a 

vHi) '^ \ - ^(7) 9 e(^) g , 



b 



(7)' 



I c 



e(7)Q<5^^ - 6(^)^(5^0, , 



(B9) 



(BIO) 



(it is perhaps worth mentioning that there is a change in 
notation, the role of Hanson and Regge's A^"^ matrix 3 
is played by the tetrad vectors [H, 0] e(^)'' here). Now, it 
is straightforward to realize that [p, $2] = 4e(o)^ W 7^ 
(if 6(0)^ W'^ = 0, the proof ends here). 
Squaring (jB6jl one gets 



(Bll) 



1 2 2 2 t2 

1 = p — rrn J . 
We may, therefore, choose the gauge 



(B12) 



which means that r = 0, thus ending the proof. 

The same result may be achieved by computing [t, <I>2] 
and choosing the gauge t ~ 0. 
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